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Abstract 

The quadratic pion scalar radius, (r 2 )^, plays an important role for present precise deter- 
minations of 7T7T scattering. Recently, Yndurain, using an Omnes representation of the null 
isospin(I) non-strange pion scalar form factor, obtains (r 2 )^ = 0.75 db 0.07 fm 2 . This value is 
larger than the one calculated by solving the corresponding Muskhelishvili-Omnes equations, 
(r 2 )J = 0.61 ± 0.04 fm 2 . A large discrepancy between both values, given the precision, then 
results. We reanalyze Yndurain's method and show that by imposing continuity of the re- 
sulting pion scalar form factor under tiny changes in the input tttt phase shifts, a zero in the 
form factor for some S-wave 1=0 T— matrices is then required. Once this is accounted for, 
the resulting value is (r 2 )^ = 0.65 ±0.05 fm 2 . The main source of error in our determination 
is present experimental uncertainties in low energy S-wave 1=0 tttt phase shifts. Another 
important contribution to our error is the not yet settled asymptotic behaviour of the phase 
of the scalar form factor from QCD. 



1 Introduction 



The scalar form factor of the pion, T n (t), corresponds to the matrix element 
Performing a Taylor expansion around t — 0, 

r.(t) = r.(o)|i + it(r 2 )i + o(t 2 )| , (i.2) 

where (r 2 )^ is the quadratic scalar radius of the pion. 

The quantity (r 2 )J contributes around 10% p] to the values of the S-wave tttt scattering lengths 
Oq and a 2 , as determined in ref.pQ, by employing Roy equations and xPT to two loops. If one 
takes into account that this reference gives a precision of 2.2% in its calculation of the scattering 
lengths, a 10% of contribution from (r 2 )^ is a large one. Related to that, (r 2 )^ is also important 
in SU (2) x SU{2) xPT since it gives the low energy constant £4 that controls the departure of F v 
from its value in the chiral limit [21 [3] at leading order correction. 

Based on one loop xPT, Gasser and Leutwyler [2] obtained (r 2 )1 = 0.55 ± 0.15 fm 2 . This 
calculation was improved later on by the same authors together with Donoghue [4], who solved 
the corresponding Muskhelishvili-Omnes equations with the coupled channels of tttt and KK. The 
update of this calculation, performed in ref.JT], gives (r 2 )1 = 0.61 ±0.04 fm 2 , where the new results 
on S-wave 1=0 tttt phase shifts from the Roy equation analysis of ref . [5] are included. Moussallam 
[6] employs the same approach and obtains values in agreement with the previous result. 

One should notice that solutions of the Muskhelishvili-Omnes equations for the scalar form 
factor rely on non-measured T— matrix elements or on assumptions about which are the channels 
that matter. Given the importance of (r 2 )^, and the possible systematic errors in the analyses 
based on Muskhelishvili-Omnes equations, other independent approaches are most welcome. In 
this respect we quote the works [3 El [9], and Yndurain's ones [TQl [TTJ [12] . These latter works have 
challenged the previous value for (r 2 )^, shifting it to the larger (r 2 )J = 0.75 ± 0.07 fm 2 . From 
ref. [TJ the equations, 

= +0.027A r2 , 5a 2 = -0.004A r 2 , (1.3) 

give the change of the scattering lengths under a variation of (r 2 )1 defined by (r 2 )1 = 0.61(1 + 
A r 2) fm 2 . For the difference between the central values of (r 2 )1 given above from refs.[H [ID] , one 
has A r 2 = +0.23. This corresponds to <5a[J = +0.006 and 5a^ = —0.001, while the errors quoted 
are = 0.220 ± 0.005 and a 2 = -0.0444 ± 0.0010. We then adduce about shifting the central 
values for the predicted scattering lengths at the level of one sigma. 

The value taken for (r 2 )J is also important for determining the 0(p A ) xPT coupling £ 4 . The 
value of ref. pQ is £4 = 4.4 + 0.2 while that of ref. [10J is £4 = 5.4 + 0.5. Both values are incompatible 
within errors. 

The papers [T01 [HI E2] have been questioned in refs.[T31 H3]- The value of the Ktt quadratic 
scalar radius, (r 2 )f n , obtained by Yndurain in ref. [ID], (r 2 )^ = 0.31 ± 0.06 fm 2 , is not accurate, 
because he relies on old experiments and on a bad parameterization of low energy S-wave 1=1/2 Ktt 
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phase shifts by assuming dominance of the k resonance as a standard Breit-Wigner pole [15]. Fur- 
thermore, (r 2 )^ was recently fixed by high statistics experiments in an interval in agreement with 
the sharp prediction of [15J, based on dispersion relations (three-channel Muskhelishvili-Omnes 
equations from the T— matrix of ref. [16]) and two-loop xPT [T7]. From the recent experiments 
[111 [[9], one has for the charged kaons [18] (r 2 )f ±7r = 0.235 ± 0.014 ± 0.007 fm 2 , and for the 
neutral ones [19] {r 2 )f L7T = 0.165 ± 0.016 fm 2 . The prediction of [15], in an isospin limit, is 
(r 2 )^ n = 0.192 ± 0.012 fm 2 , lying just in the middle of the experimental determinations. Another 
issue is Yndurain's more sound determination of the pionic scalar radius, whose (in) correctness is 
not settled yet. 

In this paper we concentrate on the approach of Yndurain [TQl HH [12] to evaluate the quadratic 
scalar radius of the pion based on an Omnes representation of the 1=0 non-strange pion scalar form 
factor. Our main conclusion will be that this approach [10] and the solution of the Muskhelishvili- 
Omnes equations [4], with irn and KK as coupled channels, agree between each other if one 
properly takes into account, for some T— matrices, the presence of a zero in the pion scalar form 
factor at energies slightly below the KK threshold. Precisely these T— matrices are those used in 
[TO] and favoured in [TTJ. Once this is considered we conclude that (r 2 )1 = 0.63 ± 0.05 fm 2 . 

The contents of the paper are organized as follows. In section 2 we discuss the Omnes rep- 
resentation of T n (t) and derive the expression to calculate (r 2 ) 1 ^. This calculation is performed 
in section 3, where we consider different parameterizations for experimental data and asymptotic 
phases for the scalar form factor. Conclusions are given in the last section. 



2 Scalar form factor 



The pion scalar form factor r„.(£), eq. fll.ll) . is an analytic function of t with a right hand cut, 
due to unitarity, for t > 4m 2 . Performing a dispersion relation of its logarithm, with the possible 
zeroes of T^t) removed, the Omnes representation results, 



r w (t) = P(t)exp 



7T 7 4m 2 S(S - t) 



ds 



(2.1) 



Here, P(t) is a polynomial made up from the zeroes of r 7r (t), with -P(O) = T^O). In the previous 
equation, 0(s) is the phase of T v (t)/ P(t), taken to be continuous and such that 0(4m 2 ) = 0. In 
ref. [10] the scalar form factor is assumed to be free of zeroes and hence P(t) is just the constant 
r„-(0) (the exponential factor is 1 for t = 0). Thus, 



r n {t) = r 7r (o)ex P 



_tt J 4m 2 s(s - t) 



ds 



From where it follows that, 
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-ds 



(2.2) 



(2.3) 



Ami 



One of the features of the pion scalar form factor of refs.jH [8], as discussed in ref. [13], is the 
presence of a strong dip at energies around the KK threshold. This feature is also shared by the 
strong S-wave 1=0 tttt amplitude, t n7T . This is so because t n7T is in very good approximation purely 
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elastic below the KK threshold and hence, neglecting inelasticity altogether in the discussion that 
follows, it is proportional to sin S n e tS7r , with 5 n the S-wave 1=0 txtt phase shift. It is an experimental 
fact that 5 n is very close to it around the KK threshold, as shown in figO Therefore, if 5 n = tt 
happens before the opening of this channel the strong amplitude has a zero at that energy. On 
the other hand, if 5 n = tt occurs after the KK threshold, because inelasticity is then substantial, 
see eq. (12.41) below, there is not a zero but a pronounced dip in \t n7T \. This dip can be arbitrarily 
close to zero if before the KK threshold 5^ approaches tt more and more, without reaching it. 
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Figure 1: S-wave I = it-it phase shift, (5 T (s). Experimental data are from refs. [2T1 [25| [26 | [27] . 

Because of Watson final state theorem the phase 0(s) in eq. ( 12. 1[) is given by 5 n (s) below the 
KK threshold, neglecting inelasticity due to 4tt or 6tt states as indicated by experiments [20J . The 
situation above the KK threshold is more involved. Let us recall that 

U, = (ve 2lS - -l)/2z , (2.4) 

with < 7] < 1 and the inelasticity is given by 1 — rj 2 , with 77 the elasticity coefficient. We denote 
by ip(s) the phase of t^, required to be continuous (below Am 2 K it is given by S n (s)). By continuity, 
close enough to the KK threshold and above it, 77 — > 1 and then we are in the same situation as in 
the elastic case. As a result, because of the Watson final state theorem and continuity, the phase 
(f>(s) must still be given by y?(s). For 5 n (sK) < tt, sk = 4?^^, ip(s) does not follow the increasing 
trend with energy of S n (s) but drops as a result of eq. (12.41) . see fig f2] for 5 n (sK') < 7r. This is easily 
seen by writing explicitely the real and imaginary parts of t n7T in eq. (l2.4p . 

1 i 

t n7T = -rism25 n H — (1 — rj cos25 n ) . (2.5) 
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Figure 2: Left panel: Strong phase <p(s), eigenvalue phase Sr + \(s) and asymptotic phase as (s). Right 
panel: Integrand of {r 2 }^ in eq. (|3.12j ) for parameterization I (dashed line) and II (solid line). For more 
details see the text. Notice that the uncertainty due to QS (s) is much reduced in the integrand. 

The imaginary part is always positive (rj < 1 above the KK threshold and 1.1 GeV [20]) while the 
real part is negative for 5 n < it, but in an interval of just a few MeV the real part turns positive 
as soon as 5 n > ir, figflj As a result, (p(s) passes quickly from values below but close to n to 
the interval [0, 7r/2]. This rapid motion of <ft(s) gives rise to a pronounced minimum of at 
this energy, as indicated in ref. [13] and shown in ngj3l The drop in 0(s) becomes more and more 
dramatic as 5 n (sx) — * (with the superscript +(— ) indicating that the limit is approached from 
values above(below), respectively); and in this limit, 0(sfc) = <p(sk) is discontinuous at sk- This is 
easily understood from eq. (12.51) . Let us call S\ the point at which 5 n (si) = n with si > sk- Close 
and above S\, (p(s) € [0, vr/2], for the reasons explained above, and <p(s) has decreased very rapidly 
from almost 7r at the KK threshold to values below 7r/2 just after s\. Then, in the limit s% — > 
one has (p{s~ bC ) = ip(s^) = n on the left, while on the right 4>{s K ) = <f(s K ) < tt/2. As a result ip(s) 
is discontinuous at s = sk- We stress that this discontinuity of ip(s) at sk when S^sk) — > k~ 
applies rigorously to 4>{sk) as well since T](sk) = 1- This discontinuity at s = sk implies also that 
the integrand in the Omnes representation for T n (t) develops a logarithmic singularity as, 

^liM! log A ; (2. 6) 

7T S K 

with 5 — > + . When exponentiating this result one has a zero for T n (sK) as (5/sk) u , v = (4>(s K ) — 
4>(s K ))/n > and 5 — > + . This zero is a necessary consequence when evolving continuously from 
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3tt(sk) < ft to 5 n (sK) > 7tCj This in turn implies rigorously that in the Omnes representation of 



T n (t), eq. (12.11) . P(t) must be a polynomial of first degree for those cases with 5 n (sK) > ft 



#2 



P(t) = 1^(0)^, (2.7) 

Si 

with si the position of the zero. Notice that the degree of the polynomial P(t) is discrete and thus 
by continuity it cannot change unless a singularity develops. This is the case when 5 n (s k) = ft, 
changing the degree from to 1. Hence, if S^^sk) > ft for a given t nn , instead of eqs. (12.21) and 
(12. 3p one must then consider, 



T n (t) = r w (0)^— ^exp 
si 



ft Jiml S(S - t) 



ds 



(2- 



and 



For those for which 6„(sk) > ft then <p(s) follows S n (s) just after the KK threshold and there 
is no drop, as emphasized in ref. [H], see figEJ 

Summarizing, we have shown that T n (t) has a zero at si when S^sk) > ft as a consequence of 
the assumption that 4>(s) follows <p(s) above the KK threshold, along the lines of ref. [TTJ, and by 
imposing continuity in T n (t) under small changes in S n (sK) — ft- As a result eqs.( 12.8p and ( 12. 9h 
should be used in the latter case, instead of eqs.( l2.2p and ( 12.31) . valid for 5 n (s k) < ft- This solution 
was overlooked in refs. [TTJ], [TH UJ] . We show in appendix A why the previous discussion on the 
zero of T v {t) for S^sk) > vr at si cannot be applied to all pion scalar form factors, in particular 
to the strange one. 

If eq.( l2.2p were used for those t n7T with ^(s^) > tt then a strong maximum of (T^t)! would 
be obtained around the KK threshold, instead of the aforementioned zero or the minimum of 
refs. |H E], as shown in figfj] by the dashed-dotted line. That is also shown in fig. 10 of ref. |22j 
or fig. 2 of [13]. This is the situation for the T n (t) of refs. [TO! [TT] . and it is the reason why (r 2 )J 
obtained there is much larger than that of refs.JU [TJ [6]. That is, Yndurain uses eqs. (l2.2p . (12. 3p for 
8n{sK) > ft, instead of eqs. (l2.8p . (12. 9p (solid line in figE]). The unique and important role played 
by ^(sx) (for elastic t nn below the KK threshold) is perfectly recognised in ref. [HJ. However, in 
this reference the astonishing conclusion that r 7r (t) has two radically different behaviours under 
tiny variations of t nn was sustained. These variations are enough to pass from 5^(3 k) < ft to 
^tt(sk) > ft [TP] , while the T— or S— matrix are fully continuous. Because of this instability of the 
solution of refs.pJJJ E] under tiny changes of 5 n (s), we consider ours, that produces continuous 
T n (t), to be certainly preferred. We also stress that our solutions, either for S n (sK) > ft and 
^{sk) < ft, are the ones that agree with those obtained by solving the Muskhelishvili-Omnes 
equations [11 [JJ E] and Unitary %PT [5] . 



can be shown from eq. (|2.5|l that 4>(s K ) — 0(s^-) = 7r. Here we are assuming rj = 1 for s < sk, which is a 
very good approximation as indicated by experiment [20\ 12 lj . 

^ 2 We are focusing in the physically relevant region of experimental allowed values for <5 t (sk), which can be larger 
or smaller than n but close to. 
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Figure 3: \T w (t) /T w (0)\ from eq. (|2.2j) with 5 1t (sk) < tt, dashed-line, and S 7t (sk) > tt, dashed-dotted line. 
The solid line corresponds to use eq. (12.8j) for the latter case. For this figure we have used parameterization 
II (defined in section [3]) with a± = 2.28 (dashed line) and 2.20 (dashed-dotted and solid lines). The 
dashed-double-dotted line is the scalar form factor of ref.jH] that has 6 w (sk) > tt- 

Let us now show how to fix si in terms of the knowledge of 5 n (s) with 5 n (sK) > 7r. For this 
purpose let us perform a dispersion relation of r 7r (t) with two subtractions, 

1 t 2 r°° ImF (?) 

r«{t) = rv(o) + Ur 2 y s t + - / =r*%d8 , (2.io) 

6 7T J Am 2 S 2 {S - t) 

From asymptotic QCD [23J one expects that the scalar form factor vanishes at infinity [101 [12] . 
then the dispersion integral in eq. fl2.10l) should converge rather fast. Eq. fl2.10p is useful because 
it tells us that the only point around 1 GeV where there can be a zero in T w (t) is at the energy 
Si for which the imaginary part of T n (t) vanishes. Otherwise, the integral in the right hand side 
of eq. (12.101) picks up an imaginary part and there is no way to cancel it as r 7r (0), (r 2 )^ and t are 
all real. Since |Imr 7r (t)| = |r^(t) sin #„.(£) | for t < s K , it certainly vanishes at the point si where 
5 n (si) = 7i. As there is only one zero at such energies, this determines s\ exactly in terms of the 
given parameterization for ^(s). 

One could argue against the argument just given to determine s\ that this energy could be 
complex. However, this would imply two zeroes at si and s*, and then the degree of P(t) would 
be two instead of one. Notice that the degree of the polynomial P(t) is discrete and thus, by 
softness in the continuous parameters of the T— matrix, its value should stay at 1 for some open 
domain in the parameters with 5k(sk) > tt until a discontinuity develops. Physically, the presence 
of two zeroes would in turn require that <p(s) — ► 3tt so as to guarantee that T n (t) still vanishes as 
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— 1/t, as required by asymptotic QCD [23], [TO] . This value for the asymptotic phase seems to be 
rather unrealistic as (f(s) only reaches 2n at already quite high energy values, as shown in figfSJ 



3 Results 



Our main result from the previous section is the sum rule to determine (r 2 ) 
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9(5n(s K ) - tt) + 




-ds 



(3.11) 



/ 2 



where 9{x) = for x < and 1 for x > 0. We split (r 2 )^ in two parts: 




Qh + Qa , 



Qh 



9(5 n (s K ) - tt) + 




ds 



77 J Ami S ^ 



Qa 



6 r° </>(*) 



(3.12) 



with s# = 2.25 GeV 2 . Reasons for fixing sh to this value are given below. 

The main issue in the application of eq. fl3.lll) is to determine <p(s) in the integrand. Below the 
KK threshold and neglecting inelasticity, one has that <j){s) = 5 T (s), 4m 2 < s < Am 2 K . This follows 
because of the Watson final state theorem, continuity and the equality 0(4m 2 ) = S n (4m 2 ) = 0. 

For practical applications we shall consider the S-wave 1=0 ittt phase shifts given by the 
K— matrix parameterization of ref.JSU] (from its energy dependent analysis of data from 0.6 GeV 
up to 1.9 GeV) and the parameterizations of ref.pQ (CGL) and ref. [21] (PY). The resulting 5 n (s) 
for all these parameterizations are shown in figO We use CGL from nir threshold up to 0.8 GeV, 
because this is the upper limit of its analysis, while PY is used up to 0.9 GeV, because at this 
energy it matches well inside the experimental errors with the data of |20j. The K— matrix of 
ref. (2D] is used for energies above 0.8 GeV, when using CGL below this energy (parameterization 
I), and above 0.9 GeV, when using PY for lower energies (parameterization II). We take the pa- 
rameterizations CGL and PY as their difference below 0.8 GeV accounts well for the experimental 
uncertainties in 5^, see figiH and they satisfy constraints from xPT (the former) and dispersion re- 
lations (both). The reason why we skip to use the parameterization of ref.|20J for lower energies is 
because one should be there as precise as possible since this region gives the largest contribution to 
(r 2 )J, as it is evident from the right panel of figfSJ It happens that the K— matrix of [20J, that fits 
data above 0.6 GeV, is not compatible with data from K ei decays [251 [26]. We show in the insert 
of figH] the comparison of the parameterizations CGL and PY with the K e ^ data of [251 126] . We 
also show in the same figure the experimental points on 5 n from refs. [20l [2T| [27]. Both refs. [201 f2Tj 
are compatible within errors, with some disagreement above 1.5 GeV. This disagreement does not 
affect our numerical results since above 1.5 GeV we do not rely on data. 

The K— matrix of ref. [2U] is given by, 



Kij(s) = ajO,-/ {x x - s) + PiPj/ (x 2 - s) + jij , 



(3.13) 



8 



where 



x{ = 0.11 ±0.15 
ai = 2.28 ±0.08 
p 1 = -1.00 ±0.03 



xl /2 = 1.19 ±0.01 
a 2 = 2.02 ±0.11 
= 0.47 ±0.05 



(3.14) 



7n = 2.86 ±0.15 7i 2 = 1.85 ± 0.18 7 22 = 1.00 ± 0.53 , 

with units given in appropriate powers of GeV. In order to calculate the contribution from the 
phase shifts of this K— matrix we generate Monte-Carlo gaussian samples, taking into account the 
errors shown in eq. fl3.14l) . and evaluate Qh according to eq. (13.12fl . The central value of 5. k (sk) 
for the K— matrix of ref. [20] is 3.05, slightly below n. When generating Monte-Carlo gaussian 
samples according to eq. fl3.14p . there are cases with S^Sk) > 7T, around 30% of the samples. Note 
that for these cases one also has the contribution — 6/sj in eq. fl3.lll) . 

The application of Watson final state theorem for s > 4m|- is not straightforward since inelastic 
channels are relevant. The first important one is the KK channel associated in turn with the 
appearance of the narrow / (980) resonance, just on top of its threshold. This implies a sudden 
drop of the elasticity parameter 77, but it again rapidly raises (the /o(980) resonance is narrow 
with a width around 30 MeV) and in the region l.l 2 < s < 1.5 2 GeV 2 is compatible within 
errors with r\ = 1 [201 EI]- F° r V — 1> the Watson final state theorem would imply again that 
4>(s) = ip(s), but, as emphasized by [13], this equality only holds, in principle, modulo tt. The 
reason advocated in ref.[13j is the presence of the region sk < s < l.l 2 GeV 2 where inelasticity 
can be large, and then continuity arguments alone cannot be applied to guarantee the equality 
~ <p(s) for s > l.l 2 GeV 2 . This argument has been proved in ref. [11] to be quite irrelevant 
in the present case. In order to show this a diagonalization of the tctc and KK S— matrix is done. 
These channels are the relevant ones when 77 is clearly different from 1, between 1 and 1.1 GeV. 
Above that energy one also has the opening of the r]f] channel and the increasing role of multipion 
states. 

We reproduce here the arguments of ref. [11], but deliver expressions directly in terms of the 
phase shifts and elasticity parameter, instead of K— matrix parameters as done in ref. [11]. For 
two channel scattering, because of unitarity, the T— matrix can be written as: 



with Sk the elastic S-wave 1=0 KK phase shift. In terms of the T-matrix the S-wave 1=0 S— matrix 
is given by, 



where the K— matrix is real and symmetric along the real axis for s > 4m 2 and Q = diag(q n , g^-), 
with q n (qK) the center of mass momentum of pions(kaons). This allows one to diagonalize K with 
a real orthogonal matrix C, and hence both the T— and S— matrices are also diagonalized with 
the same matrix. Writing, 




(3.15) 



S = I + 2iT , 

satisfying SS* = S^S = I. The T-matrix can also be written as 

T = Q 1 / 2 (K- 1 - tQ)- 1 Q 1 ' 2 , 



(3.16) 



(3.17) 




(3.18) 
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one has 



cos 6 



sin 9 



/2 



1/2 



1 — r] 2 cos 2 A — 77 1 sin A | \Jl — rf cos 2 A 
sin A r] - + (1 - rf) cot 2 A 

v 1 — rf cos 2 A — sin A| yl — r? 2 cos 2 A 

with A = 5k — S n . On the other hand, the eigenvalues of the S— matrix are given by, 



(3.19) 



e 2i5 (+) 



e 2j5 (-) 



Si 



1 + e 



2iA 



1 - - tan A \jl + (1 - r] 2 ) cot 2 A 



5 



1 + e 



-2iA 



22- 



1 + - tan A J 1 + (1 - T] 2 ) cot 2 A 



(3.20) 
(3.21) 



The eigenvalue phase <5( + ) satisfies 5( + )(sx) = 8tt( s k)- The expressions above for exp2z5( + ) and 
exp 2z5(_) interchange between each other when tan A crosses zero and simultaneously the sign in 
the right hand side of eq. fl3.19l) for sin# changes. This diagonalization allows to disentangle two 
elastic scattering channels. The scalar form factors attached to every of these channels, 1^ and 
r' 2 , will satisfy the Watson final state theorem in the whole energy range and then one has, 



r' 



= ( rl) = ° Tq1/2t = ° TQl/2 ( rv ) ' 

= q- 1/2 (Acos^|r / 1 |e l<5 (+) ±sm6\T' 2 \e iS (-->) , 
= q K 1/2 (±cos6\T' 2 \e l5 (--> - A sin 9 |ri|e w c+>) 



(3.22) 



The ± in front of \T' 2 \ is due to the fact that T' 2 



at Sk, as follows from its definition in the 
equation above. Since Watson final state theorem only fixes the phase of T 2 up to modulo tt, and 
the phase is not defined in the zero, we cannot fix the sign in front at this stage. Next, 1^ has a 
zero at s x when ^(s^) > tt. For this case, —\T[\ must appear in the previous equation, so as to 
guarantee continuity of its ascribed phase, and this is why A = (— l) e ( <5 ' r (' s ^) _7r ). 

Now, when rj — > 1 then sin^ — > as — T])/2 and (J)(s) is then the eigenvalue phase <5( + ). 
This eigenvalue phase can be calculated given the T— matrix. For those T— matrices employed 
here, and those of refs. [TU l fTT l HI [T3] . <S(+)(s) follows rather closely (p(s) in the whole energy range. 
This is shown in figj2] and already discussed in detail in ref . [11] . In this way, one guarantees 
that <f)(s) and (p(s) do not differ between each other in an integer multiple of tt when 77 ~ 1, 
l.l 2 < 8 < 1.5 2 GeV 2 . 

For the calculation of Qh in eq. (I3.12j) we shall equate 0(s) = ip(s) for Amj^ < s < 1.5 2 GeV 2 . 
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Denoting, 



iff 
h 

h 
h 



6 



^Y = h + h + h, 



7T ,/4 m 2 S 

6 r K 



7T ,/4 m 2 S 



2 



l.V 



c/.s 



then 



6 

— 6>(^(sk) 
si 



7T 



(3.23) 



(3.24) 



Now, eq. fl3.22p can also be used to estimate the error of approximating 0(s) by ip(s) in the range 
Am 2 K < s < 1.5 2 GeV 2 to calculate 7 2 and 7 3 as done in eq. (13.23 p . We could have also used 5(+)(s) 
in eq. (l3.23p . However, notice that when 77 < 1 then <p(s) ~ and when inelasticity could be 

substantial the difference between <5( + )(s) and <p(s) is well taken into account in the error analysis 
that follows. Remarkably, consistency of our approach also requires <f>(s) to be closer to <p(s) than 
to 8(+\(s). The reason is that <p(s) for ^(s^) > n is in very good approximation the <p(s) for 
£>k{sk) < K plus ^1 this is clear from figfSJ This difference is precisely the required one in order 
to have the same value for (r 2 )1 either for 5^(sk) < 7r or ^(s^) > n from eq. (13. llj) . However, the 
difference for 5(+)(s) between S^^sk) < vr and S w (sk) > 7r is smaller than ir. Indeed, we note that 
4>(s) follows closer ip(s) than <5(+)(s) for the explicit form factors of refs.[8j|I]. 

Let us consider first the range l.l 2 < s < 1.5 2 GeV 2 where from experiment [2U] r\ — 1 within 



errors. With e = ± tan ^Ir^/r'jj and p 



*<-) 



eq. (13.221) allows us to write, 



I\ = Acos0|ri|e w (+)(l + ecosp) ( 1 + i 



e sm p 
1 + e cos p 



(3.25) 



When 77 — > 1 then e — > 0, according to the expansion 



#3 



tan# 



sin A 



' l + 3cos2A, i 

1 a • 2 a (i-v) 

8 sm A 



Rewriting, 



e sm p 



exp 1 



e sm p 



0(e 2 ) 



1 + ecosp \ l + ecosp / 

which from eqs. (13.251) and (13.271) implies a shift in <5( + ) because of inelasticity effects, 



'(+) °(+) 



Sm + 



e sm p 
1 + e cos p 



(3.26) 



(3.27) 



(3.28) 



# 3 The the ratio ir^/r^J, present in e, is not expected to be large since the /o(1300) couples mostly to 4ir and 
similarly to irir and KK, and the /o(1500) does mostly to mr [25] , 
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Using rj = 0.8 in the range l.l 2 < s < 1.5 2 GeV, 77 ~ 1 from the energy dependent analysis of 
ref. [2n] given by the K— matrix of eq.f l3.13l) . one ends with e ~ 0.3. Taking into account that <5( + ) 
is larger than > 3n/2 for 8 w (sk) > tt (in this case ~ (5^) , and around 37r/4 for 5 n (sK) < tt, 
see figJ2J one ends with relative corrections to <5( + ) around 6% for the former case and 13% for the 
latter. Although the K— matrix of ref. [20], eq. (l3.13[) . is given up to 1.9 GeV, one should be aware 
that to take only the two channels tttt and KK in the whole energy range is an oversimplification, 
particularly above 1.2 GeV. Because of this we finally double the previous estimate. Hence ^3 is 
calculated with a relative error of 12% for 5 7T (sk) > tt and 25% for 8 n (spc) < tt. 

In the narrow region between sk < s < l.l 2 GeV 2 , rj can be rather different from 1, due to 
the /o(980) that couples very strongly to the just open KK channel. However, from the direct 



measurements of ttti — ► KK [29], where 1 — rj 2 is directly measuredP 4 ! one has a better way to 
determine rj than from txtt scattering [201 121] . It results from the former experiments, as shown 
also by explicit calculations J3UIEI1E2], that rj is not so small as indicated in im experiments [2"U] . 
and one has rj ~ 0.6 — 0.7 for its minimum value. Employing rj = 0.6 in eq. (13.281) then e ~ 0.5. 
Taking around 7r/2 when S^(sk) < n this implies a relative error of 30%. For 8^{sk) > it one 
has instead 8<+\ > it, and a 15% of estimated error. Regarding the ratio of the moduli of form 
factors entering in e we expect it to be < 1 (see appendix A). Therefore, our error in the evaluation 
of I2 is estimated to be 30% and 15% for the cases S^sk) < tt and S^sk) > 7T, respectively. 

As a result of the discussion following eq. (13.241) . we consider that the error estimates done 
for I 2 and ^3 in the case 5 n (sK) < tt are too conservative and that the relative errors given for 
5tt(sk) > K ar e more realistic. Nonetheless, since the absolute errors that one obtains for I 2 and 
^3 are the same in both cases (because I2 and ^3 for 5 n {sx) < ti are around a factor 2 smaller than 
those for 5 7T (sk) > 7r) we keep the errors as given above. To the previous errors for I 2 and ^3 due 
to inelasticity, we also add in quadrature the noise in the calculation of Qh due to the error in 
t n7T from the uncertainties in the parameters of the K— matrix eqs. (l3.13p . (I3.14p . and those in the 
parameterizations CGL and PY. 

We finally employ for s > 2.25 GeV 2 the knowledge of the asymptotic phase of the pion scalar 
form factor in order to evaluate Qa in eq. (13. 121) . The function 0(s) is determined so as to match 
with the asymptotic behaviour of T n (t) as —1/t from QCD. The Omnes representation of the 
scalar form factor, eqs. (12.21) and (12. 8p . tends to t~ q ^ and t~ q ^ +l for t — > 00, respectively. Here, q 
is the asymptotic value of the phase <f>(s) when s — » 00. Hence, for S n (sj() < n the function <f>(s) 
is then required to tend to n while for 5 w (sk) > tt the asymptotic value should be 2n. The way 
(p(s) is predicted to approach the limiting value is somewhat ambiguous |TT|, [T2]. 

^ s ^*{ ii± t£w>) ■ (3 ' 29) 

In this equation, 2d m = 24/(33 — 2nf) ~ 1, A 2 is the QCD scale parameter and n = 1, 2 for 
8 n {4m 2 K ) < it, > ii, respectively. The case n = 2 was not discussed in refs. [TO] [TT] [T2l [13] [T4] for 
the form factor given in eq. fll.ip . There is as well a controversy between [T4] and [12] regarding 
the ± sign in eq. (13.291) . If leading twist contributions dominate [TU [12] then the limiting value is 
reached from above and one has the plus sign, while if twist three contributions are the dominant 
ones p3j the minus sign has to be considered [12]. In the left panel of figf5]we show with the wide 



# 4 Neglecting multipion states. 
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1 / \ 


I 


I 


TT 

II 


TT 

II 


On{S K ) 


> TT 


< TT 


> TT 


< TT 


h 


0.435 ± 0.013 


0.435 ± 0.013 


0.483 ±0.013 


0.483 ±0.013 


h 


0.063 ±0.010 


0.020 ±0.006 


0.063 ±0.010 


0.020 ±0.006 


h 


0.143 ±0.017 


0.053 ±0.013 


0.143 ±0.017 


0.053 ±0.013 


Qh 


0.403 ± 0.024 


0.508 ±0.019 


0.452 ± 0.024 


0.554 ±0.019 


Qa 


0.21 ±0.03 


0.10 ±0.03 


0.21 ±0.03 


0.10 ±0.03 


(r 2 y s 


0.61 ±0.04 


0.61 ±0.04 


0.66 ±0.04 


0.66 ±0.04 



Table 1: Different contributions to (r 2 )^ as defined in eqs.f l3.12"l) and (13.231) . All the units are fm 2 . 
In the value for {r 2 )1 the errors due to Ji, I2, h and Qa are added in quadrature. 



bands the values of (p(s) as for s > 2.25 GeV 2 from eq. fl3.29p . considering both signs, for n — 1 
(S^sk) < tt) and 2 (^( Sjr ) > vr). We see in the figure that above 1.4-1.5 GeV (1.96-2.25 GeV 2 ) 
both (p(s) and 4>(s) as phases match and this is why we take sh = 2.25 GeV 2 in eq.f l3.lll) . similarly 
as done in refs.|10[ [TT]. In this way, we also avoid to enter into hadronic details in a region where 
n < 1 with the onset of the /o(1500) resonance. The present uncertainty whether the + or — sign 
holds in eq. (13.291) is taken as a source of error in evaluating Qa- The other source of uncertainty 
comes from the value taken for A 2 , 0.1 < A 2 < 0.35 GeV 2 , as suggested in ref.[lOJ. From figj2]it is 
clear that our error estimate for <p as (s) is very conservative and should account for uncertainties 
due to the onset of inelasticity for energies above 1.4 — 1.5 GeV and to the appearance of the 
/o(1500) resonance. In the right panel of figj2] we show the integrand for (r 2 )^, eq. fl3.12p . for 
parameterization I (dashed line) and II (solid line). Notice as the large uncertainty in <p a s{s) is 
much reduced in the integrand as it happens for the higher energy domain. 

In tabled] we show the values of Ji, I2, h, Qh, Qa and (r 2 )J for the parameterizations I and 
II and for the two cases S n (s k) > n and S^sk) < tt. This table shows the disappearance of the 
disagreement between the cases S^sk) > tt and S n (sK) < tt from the tttt and KK T— matrix 
of eq. fl3.13p . once the zero of T^t) at si < sk is taken into account for the former case. This 
disagreement was the reason for the controversy between Yndurain and ref . [13] regarding the value 
of (r 2 )1 . The fact that the parameterization II gives rise to a larger value of (r 2 )1 than I is because 
PY follows the upper 6^ data below 0.9 GeV, while CGL follows lower ones, as shown in figfTJ 

The different errors in table [1] are added in quadrature. The final value for (r 2 )J is the mean 
between those of parameterizations I and II and the error is taken such that it spans the interval 
of values in table [1] at the level of two sigmas. One ends with: 

(r 2 y s = 0.63 ± 0.05 fm 2 . (3.30) 

The largest sources of error in (r 2 )1 are the uncertainties in the experimental 5 n and in the 
asymptotic phase <p as . This is due to the fact that the former are enhanced because of its weight 
in the integrand, see figJ21 and the latter due to its large size. 

Our number above and that of refs. [JJ, H] , (r 2 )1 = 0.61 ±0.04 fm 2 , are then compatible. On the 
other hand, we have also evaluated (r 2 )1 directly from the scalar form factor obtained with the 
dynamical approach of ref. [8] from Unitary xPT and we obtain (r 2 )J = 0.64 ±0.06 fm 2 , in perfect 
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agreement with eq. fl3.30p . Notice that the scalar form factor of ref. [8] has S w (sk) > n and we have 
checked that it has a zero at si, as it should. This is shown in figj3]by the dashed-double-dotted line. 
The value (r 2 )1 = 0.75 ±0.07 fm 2 from refs. |10[ [TT] is much larger than ours because the possibility 
of a zero at Si was not taking into account there and other solution was considered. This solution, 
however, has an unstable behaviour under the transition ^(s^) = 7r — + to ^(s^) = 7r + + and 
it cannot be connected continuously with the one for 5 7T (sk) < 7r. Our solution for T n (t) from 
Yndurain's method does not have this unstable behaviour and it is continuous under changes in 
the values of the parameters of the K— matrix, eqs. fl3.13p and (13.141) . This is why, from our results, 
it follows too that the interesting discussion of ref.pj], regarding whether 5 n (sK) < n or > ir, 
is not any longer conclusive to explain the disagreement between the values of refs. [TOT [TT] and 
ref.pQ for (r 2 ) w s . 

We can also work out from our determination of (r 2 )^, eq.f l3.30l) . values for the 0(p 4 ) SU{2) 
xPT low energy constant £4. We take the two loop expression in xPT for (r 2 )J [TJ, 

where f n = 92.4 MeV is the pion decay constant, £ = (M n / 4tt f v ) 2 and M n is the pion mass. First, 
at the one loop level calculation A r = and then one obtains, 

h = 4.7 ±0.3 . (3.32) 

We now move to the determination of £4 based on the full two loop relation between (r 2 )1 and £4. 
The expression for A r can be found in Appendix C of ref.pQ. A r is given in terms of one 0(p 6 ) 
xPT counterterm, r S2 , and four 0(p 4 ) ones. Taking the values of all these parameters, but for £4, 
from ref.|T], and solving for £4, one arrives to 

h = 4.5 ± 0.3 . (3.33) 

This number is in good agreement with £4 = 4.4 ± 0.2 [TJ. 

Ref.[T2] also points out that one loop xPT fits to the S-, P- and D-wave scattering lengths and 
effective ranges give rise to much larger values for £ 2 and £4 than those of ref . [JJ . For more details 
we refer to [T21 . 



4 Conclusions 

In this paper we have addressed the issue of the discrepancies between the values of the quadratic 
pion scalar radius of Leutwyler et al. jU [13], {r 2 )1 = 0.61 ± 0.04 fm 2 , and Yndurain's papers 
PHI CHI H2] , (r 2 ) = 0.75±0.07 fm 2 . One of the reasons of interest for having a precise determination 
of (r 2 )g is its contribution of a 10% to a[} and a 2 ,, calculated with a precision of 2% in ref. [TJ. The 
value taken for (r 2 )1 is also important for determining the 0{p A ) %PT coupling £ 4 . 

From our study it follows that Yndurain's method to calculate (r 2 )^ [101 [UJ, based on an 
Omnes representation of the pion scalar form factor, and that derived by solving the two(three) 
coupled channel Muskhelishvili- Omnes equations [H [TJ [6], are compatible. It is shown that the 
reason for the aforementioned discrepancy is the presence of a zero in T n (t) for those S-wave 1=0 
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T— matrices with S n (sx) > tt and elastic below the KK threshold, with sk = ^m? K . This zero 
was overlooked in refs. [TUt [TT] . though, if one imposes continuity in the solution obtained under 
tiny changes of the tttt phase shifts employed, it is necessarily required by the approach followed 
there. Once this zero is taken into account the same value for (r 2 )^ is obtained irrespectively 
of whether 5 w (sk) > 7r or S^(sk) < it. Our final result is (r 2 )J = 0.63 ± 0.05 fm 2 . The error 
estimated takes into account experimental uncertainty in the values of S n (s), inelasticity effects 
and present ignorance in the way the phase of the form factor approaches its asymptotic value it, 
as predicted from QCD. Employing our value for (r 2 )J we calculate £4 = 4.5 ± 0.3. The values 
(r 2 )J = 0.61 ± 0.04 fm 2 and £4 = 4.5 ± 0.3 of ref.pQ are then in good agreement with ours. 
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Appendices 

A Coupled channel dynamics 

We take tttt and KK coupled channels and denote by F\ and F2 their respective 1=0 scalar form 
factors. Unitarity requires, 

2 

ImF, = ^F^(t- S ;)^, (A.l) 
i=i 

where \\tij\\ is the 1=0 S-wave T— matrix, is the threshold energy square of channel i and 
Pi = Qi/87T^/s, with qi its center of mass three momentum. 
A general solution to the previous equations is given by, 

F = TG ,F=(£) ,C=(g) , (A.2) 

where the functions Gi(t) do not have right hand cut. This equation is interesting as tells us that 
if pion dynamics dominate, \G\\ >> \G2\1 then F% ~ G\t\\ and the form factor phase 0(s) follows 
(p(s). As a result, like in, it has a zero at s\ below the KK threshold for 5 n (s}() > 7r, as shown 
in section [31 On the other hand, if kaon dynamics dominates, [G^^j >> \G\\, then F\ ~ G2ii2 and 
<p(s) follows the phase of ii 2 , that above the KK threshold is clearly above it. This is why for 
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the pion strange scalar form factor there is no zero at Si < sk for 5 n (s K ) > tc, indeed there is a 
maximum like that shown in figj3]by the dashed-dotted line. 

As in section [3] we now proceed to the diagnolization above the KK threshold of the renormal- 
ized T— matrix T' 



r-W. ,-($). f-*rc- ft I) 



t u = sin5 (+) e i(5 (+) , t 22 = sm6 { - } e l5 ^ . (A.3) 
The corresponding diagonal form factors F[ and F 2 , collected in the vector F', are 

~ , / i cos 9 pi 1 ^ 2 Gi — sin 9 p 9 1 ^ 2 G 2 Wn \ 

F = cV' 2 f = rcV" 2 G= . /' ' S V - " • (A.4) 

\ <sm^p 1 Gi + cos#p 2 G 2 >t 22 / 

The previous expressions allow to obtaining Fi directly in terms of the eigenphases and with clean 
separation between pion, G\, and kaon dynamics, G 2 . From eq. (13.221) it follows that, 

Fi = I cos 2 9 p^Gi — cos 9 sin 9 p 2 X ^ P\ 1 G 2 \ tu 

-1/2-1/2 1- (A ' 5) 

+ < sin 2 9 p^Gi + cos 9 sin 9 p 2 ! p 1 ^ 2 G2 [ ^22 • 

For ^(sa') > vr typical values, somewhat above the KK threshold, are e 24<5 (+> ~ +i, e 2l5{ -'> ~ — z 
and sin# > 0. For dominance of G\ one has F1/G1 ~ p^ 1 ^ + cos2#)/2 while for dominance 

]/2 1 /2 

of G2 the result is F1/G2 — — sin6> cos9p 2 p 1 < 0. The factors Gi,2 do not introduce 
any change in <p(s) with respect to its value before the opening of the KK threshold since they 



are smooth functions in scj In both cases the phase <f)(s) is larger than tc and Fi follows the 
upper trend of phases shown in figj2] (note that in this case tu is in the first quadrant though 
5-n > 7r). Now, doing the same exercise for 5 7T {sk) < vr, one has the typical values e 2w5 (+> ~ —i, 
e 2«<5 ( _) ^ _|_^ an( j s j n < q p or pj ori dominance then F\/Gi ~ Pi 1 ^ — cos2#)/2 and for the kaon 

j /2 1/2 

one F1/G2 — + sin9cos9p 2 p x < 0. Thus, in the former case the phase is > 7T/2, and follows 
the lower trend of phases of figl21 while in the latter is > tt and follows again the upper trend (this 
is the case of the strange scalar form factor) . 

The demonstration in section [3] that (^{sk) is discontinuous in the limit S n (sK) — * tt - by taking 
s\ — > s^, cannot be applied in the case of kaon dominance (e.g. pion strange scalar form factor). 
From eq. ( 1A.5I) it follows that, 

F 1 (t)~ -cos 9 sin 9 p 2 1/2 Pl 1/2 G 2 {tn- t 22 ) ■ (A.6) 

The point is that £22 for t > S\ (si — > s^) is of size comparable with that of tu (both tend to zero) 
and the phase does not follow This is not the case for pion dominance because for s± —>■ 
then sin 2 ^ — > 0, F 1 (t) ~ cos 2 9 p^Gjn, eq. (IA.5l) . and 0(s) follows 

From eq. (1A.4p we can also write |r 2 /r 2 | ~ \tu tan#/t 22 | for the case of pion dominance. Since 
typically \tu/t 22 \ — 1, as shown above for energies somewhat above the KK threshold, then 
ir'rj/r^l ~ I tan 6*| < 1. This is why we consider that equating it to 1 in section [3] is a conservative 
estimate. 



# 5 Due to the Adler zeroes this is not necessarily case close to the irir threshold. 
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